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We evaluate the suppression of J/ip production in an equi- 
librating quark gluon plasma for two competing mechanisms: 
Debye screening of colour interaction and dissociation due to 
energetic gluons. Results are obtained for S + S and Au + Au 
collisions at RHIC and LHC energies. At RHIC energies the 
gluonic dissociation of the charmonium is found to be equally 
important for both the systems while the screening of the in- 
teraction plays a significant role only for the larger systems. 
At LHC energies the Debye mechanism is found to dominate 
for both the systems. While considering the suppression of 
directly produced T at LHC energies, we find that only the 
gluonic dissociation mechanism comes into play for the ini- 
tial conditions taken from the self screened parton cascade 
model in these studies. Thus we find that a systematic study 
of quarkonium suppression for systems of varying dimensions 
can help identify the source and the extent of the suppression. 

PACS numbers: 12.38M 



Relativistic heavy ion collision experiments at the 
CERN SPS are believed W to have led to a production of 
quark gluon plasma - which existed in the early universe 
and which may be present in the core of neutron stars. 
The last two decades have seen a hectic activity towards 
identifying unique signatures of the quark-hadron phase 
transition. The suppression of J/tp production in such 
collisions has been one of the most hotly debated signals 
in this connection. 

The heavy quark pair leading to the J/ip mesons are 
produced in such collisions on a very short time-scale 
~ l/2m c , where m c is the mass of the charm quark. The 
pair develops into the physical resonance over a forma- 
tion time 7^> and traverses the plasma and (later) the 
hadronic matter before leaving the interacting system to 
decay (into a dimuon) to be detected. This long 'trek' 
inside the interacting system is fairly 'hazardous' for the 
J/ip. Even before the resonance is formed it may be 
absorbed by the nucleons streaming past it 0]. By the 
time the resonance is formed, the screening of the colour 
forces in the plasma may be sufficient to inhibit a binding 
of the cc HH . Or an energetic gluon || or a comoving 
hadron could dissociate the resonance(s). The extent 
of absorption will be decided by a competition between 
the momentum of the J/ip and the rate of expansion and 
cooling of the plasma, making it sensitive to such details 
as the speed of sound J7||8) . Thus a study of J/ip produc- 
tion is poised to provide a wealth of information about 
the evolution of the plasma and its properties. 

It has been shown M that the nucleonic absorp- 
tion (the "normal absorption"), operating on the pre- 



resonance- which is yet to evolve into a physical particle- 
is identical for J/ip, ip' , and \c- This absorption is al- 
ways present and is brought about by the nucleons (or the 
Lorentz-contracted partonic clouds) streaming past the 
pre-resonances, as mentioned earlier. A reliable quanti- 
tative estimate within Glauber model is available Q for 
this. 

In the present work we concentrate on the dissociation 
of the charmonium in quark gluon plasma due to colour 
screening and scattering with gluons and ask whether we 
can distinguish between the two mechanisms. We empha- 
size that these mechanisms are in addition to nucleonic 
absorption mentioned earlier. 

In principle the colour screening is a collective effect, 
where the presence of a large number of colour quanta 
modifies the force between c and c so that, above the 
critical temperature (T c ~ 200 MeV), we have: 

V(r) = — a/r + ar — > V(r) = — cvcxp(— ^ir>r)/r (1) 

where a and a (the string tension) are phenomenological 
parameters and /i£> is the Debye mass. 

Thus, e.g., the direct production of the J/ip is inhib- 
ited once the Debye mass is more than 0.7 GeV Eel. The 
gluonic dissociation, on the other hand, is always possi- 
ble as long as an energetic gluon can be found. They 
can always be present in the tail of the thermal distribu- 
tions and thus given sufficient time, a J/ip can always be 
dissociated in a plasma of any temperature! 

Of-course in actual practice the QGP will expand and 
cool and undergo hadronization below the critical tem- 
perature T c , and thus the hot medium will have only a 
finite life-time. This enriches the competition between 
the mechanisms of the gluonic dissociation and the De- 
bye screening for the charmonium suppression. In the 
present work we show that this also provides us with a 
handle to decipher the extent to which each mechanism 
contributes to the suppression of J /ip. 

Let us assume that a thermally equilibrated plasma is 
formed in relativistic heavy ion collisions at some time Ti 
and that the elastic scattering among the partons is suf- 
ficiently rapid to maintain thermal equilibrium. A large 
number of studies [[11],[L2| have indicated that the plasma 
thus produced may not be in a state of chemical equi- 
librium and that the quark and gluon fugacities are less 
than unity. We assume that the chemical equilibration 
proceeds dominantly via 

99 <-*■ 999, 99 <-* M- ( 2 ) 

Assuming the evolution to proceed according to Bjorken 
hydrodynamics, the evolution of the parton densities are 
given by pi: 
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where Oi = 16C(3)/tt 2 w 1.95, a 2 = 8tt 2 /15 « 5.26, 
&i = OCCS)^//^ 2 w 2.20, and b 2 = 7TT 2 N f /20 « 6.9. The 
expressions for the density and velocity weighted reaction 
rates, 
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can be found in Ref . |l3J . 

The results for the time evolution of the fugacities 
and the temperature for the initial conditions obtained 
from the self screened parton cascade model [O for 
Au + Au collisions at RHIC and LHC energies are given 
in Ref. p3|. For the S + S collisions we assume that while 
the initial fugacities are same as those for the Au + Au 
system, the initial temperatures are estimated by assum- 
ing that it scales as Tj ~ A" ". This is motivated by a 
recent study on the basis of parton saturation |15] which 
also suggests that the initial number density divided by 
Tf is nearly independent of the mass-number of the nu- 
clei. This, we believe, provides a useful initial guess, 
even though the conditions envisaged for self screening 
are not strictly met for S + S at RHIC. For the sake 
of completeness, we have given the initial conditions in 
Table 1. It may be noted that these are different from 
those used in Ref. ||, which were 'inspired' by the HI- 
JING model and which had, for example, much smaller 
fugacities. (We have verified that our computer program 
fully reproduced the results of Ref. IH], with the initial 
conditions given there.) 

We shall also introduce a energy density profile such 
that, 

e(n,r) = (l + /3)<e«>(l- r 2 /R 2 f®(R - r) (7) 

where (3 = 1/2, R is the transverse dimension of the sys- 
tem and r is the transverse distance, and < ej > is the 
energy density obtained by taking the initial tempera- 



n. The profile plays an 
boundary of the hot and 



ture as Xi and fugacities as A; 
important role in defining the 
dense dcconfincd matter. 

Having obtained the density of the partons we estimate 
the Debye mass of the medium as 



A> 



k 2 x 47ra s (A s + N f X q /6) T 2 
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where we have arbitrarily taken k as 1.5 to account for 
the corrections [ttql to the lowest order perturbative QCD 



which provides the above expression for k = 1. Results 
for other values of k are easily obtained. We shall assume 
that the J/ip can not be formed in the region where /j,d 
is more than 0.7 GeV. We can then estimate the survival 
probability of the directly produced J/ip as a function 
of its transverse momentum px by proceeding along the 
lines of Ref. @,|0. 

In order to estimate the gluonic dissociation we re- 
call |17j that the short range properties of the QCD can 
be used to derive the gluon- J/ip cross-section as: 
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where q° is the gluon energy in the rest-frame of J/ip 
and eo is the binding energy of the J ftp. The expression 
for the thermal average of this cross-section < v Te \cr > 
is given in Ref. |J. (See, also Ref. jig] for an interesting 
alternative approach.) 

We wish to have a quantitative comparison of these 
two processes and therefore it is imperative that we com- 
pare their results for similar conditions. Thus, exactly 
as while dealing with Debye screening, we assume that 
the cc produced intially takes a finite amount of time 
~ 0.89 fm/c in its rest frame to evolve into the physi- 
cal resonance. This can get large due to time dilation, 
in the frame of the plasma, leading to the characteristic 
Pt dependence of the survival probability for the J/ip 
discussed in the literature. 

We argue that the gluon- J/ip cross-section also attains 
its full value only after the cc pair has evolved into the 
physical resonance. We assume that this evolution of the 
cross-section can be parametrized as 



Co {t/t^Y if T < T,p 
CT if T > T^ 



(10) 



similarly to the case when the nuclear absorption is con- 
sidered p9| , where <jq is the cross-section estimated ear- 
lier (E q.pp T A similar assumption was invoked by Farrar 
et al. p0|[ when the QQ-system evolves as it moves away 
from the point of hard interaction. One may imagine that 
this amounts to assuming that the effective cross-section 
scales as the transverse area of the system relative to 
the size it attains when it is fully formed. In the present 
work we follow, Blaizot and Ollitrault [[19| who have used 
v = 2. Farrar et al. plj have suggested that v = 1 corre- 
sponds to a quantum diffusion of the quarks while v = 2 
would correspond to maximal rapid (classical) expansion. 
Legrand et al. |22J have used v = 1 in a recent study. 

This aspect is in contrast to the work of Xu et al || 
where a fully formed J/tp is assumed to exist right at the 
initial time in the plasma. We shall see that ignoring the 
formation time leads to an enhanced suppression of the 
charmonium. 

We can now easily estimate the time spent by the J/ip 
in the deconfined medium for a given px and get the 
survival probability following Ref. H . 



In Fig. 1 we show our results for RHIC energies for 
S + S and Au + Au collisions. We see that the combina- 
tion of a finite formation time and (reasonably) large fijj 
required to inhibit the formation of the directly produced 
J /if) in the plasma ensures that the mechanism of Debye 
screening is not effective in suppressing its production. 
However the gluonic dissociation leads to a suppression 
of the J/tp formation even after the moderating effect of 
the inclusion of formation time is included. 

The situation for the larger (and hotter) volume of 
plasma produced in Au + Au collisions is much richer in 
detail. We see that while the J / -0s having lower trans- 
verse momenta are more strongly suppressed due to the 
Debye mechanism, those having higher transverse mo- 
menta are more suppressed by the mechanism of gluonic 
dissociation. In fact we see that while the Debye screen- 
ing has become quite ineffective for p? >6 GeV, the glu- 
onic dissociation continues to be operative. The differ- 
ent results obtained here compared to authors of Ref. ||] 
(when the formation time considerations are ignored) are 
solely due to the SSPC initial conditions (Table 1) used 
here. 

The corresponding results at LHC energies are shown 
in Fig. 2. Now we see that the Debye screening is more 
effective in suppressing the production of the directly pro- 
duced J/tp at all the momenta considered, provided we 
include the considerations of the formation time while 
evaluating the gluonic dissociation, for both the systems. 

Of-course in a model calculation one can arbitrarily 
enhance the impact of Debye screening by taking a larger 
value for the coefficient « (Eq.ph. This sensitivity would 
be useful for determining its precise value J23J. 

The treatment outlined here can be extended to the 
case of T production studied in great detail by the au- 
thors of Ref. Jq,p3[, for example. We give the results 
only for the LHC energies, for the directly produced T 
(Fig. 3). We find that both for the light as well as the 
heavy systems the Debye mechanism is not at all able to 
inhibit the formation of the directly produced Ts, though 
the gluonic dissociation leads to a considerable suppres- 
sion, with the changes brought about by the inclusion 
of the formation time seen earlier for the J/ips. This 
is easily seen to be the consequence of the initial condi- 
tions used here, which have chemically non-equilibrated 
plasma leading to small Debye mass, even though the 
temperature is rather high. By the time the T is formed 
the Debye mass drops below the value of ~ 1.6 GeV, 
required to inhibit its formation, causing it to escape un- 
scathed. 

Before summarizing, let us discuss some of the assump- 
tions made in this work. We have, so far considered only 
the dissociation of the directly produced J/tp. Of-course, 
it is well known that up to 30% of the J/tps seen in these 
studies may be produced from a decay of Xc and up to 
10% or so may come from the decay of tp' , which however 
is also easily dissociated by a moderately hot (confined) 
hadronic matter and is unable to play a decisive role in 
distinguishing confined matter from a deconfined mat- 



ter. In order to include the effect of these resonances, 
we should also have access to g — Xc cross-section. This 
would involve extending the method of Ref. |p3 to the 
case of charmonium in \p state. However this is not quite 
easy as the assumption eo 3> Aqqq used in the above 
reference are not strictly valid for this case, as the binding 
energy of \c is only about 240 MeV. 

Still, from the considerations of binding energy alone, 
one may expect the g — Xc cross-sections to be larger. 
However, the Debye mass required to inhibit the forma- 
tion of Xc is also smaller and thus this competition be- 
tween the two mechanisms will continue. The inclusion 
of the transverse expansion and the possibility of a differ- 
ent speed of sound than the value of l/v3 assumed here 
will also add to the richness of the information likely to 
be available from such studies. Of-course a full study 
will additionally include the effect of the nuclear and the 
co-mover absorption, before these interesting details are 
investigated. 

The incorporation of the formation time is interesting 
for one more reason. The pre-equilibrium stage (before 
the time Tj) may be marked by presence of gluons of 
high transverse momenta, as a result of first hard colli- 
sions, and one may imagine that they play an important 
role in suppression of charmonium formation. This is 
unlikely for two reasons. Firstly, the gluon-charmonium 
cross-section drops rapidly as the gluon momentum in- 
creases 1 24 after reaching a peak around p ~ 1 GeV. 
Secondly we expect these cross-sections to be further sup- 
pressed during the formation era due to the considera- 
tions of the formation time. 

While considering the suppression of T, we found that 
only the mechanism of gluonic dissociation is playing a 
role. This happens as the initial conditions used here in- 
volve a chemically non-equilibrated plasma. If the initial 
fugacities were to be larger, the Debye screening would 
also play a role, which will definitely be a good check on 
these. 

In brief, we have seen that while the gluonic dissocia- 
tion of the J/tp is always possible, the Debye screening 
is not effective in the case of small systems at RHIC en- 
ergies. For the larger systems, the Debye screening is 
more effective for lower transverse momenta, while the 
gluonic dissociation dominates for larger transverse mo- 
menta. At LHC energies the Debye screening is the dom- 
inant mechanism of J/tp suppression for all the cases and 
momenta studied. We have also seen that the inclusion of 
the formation time of the J/tp plays an interesting role 
in reducing the role of the gluonic dissociation. As an 
interesting result, we find the gluonic dissociation to be 
substantial but the Debye screening to be ineffective for 
T suppression at the LHC energy. This may of-course 
change if different initial conditions and screening crite- 
ria P,E3[ are employed. 
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TABLE I. Initial values for the time, temperature, fugaci- 
ties etc. for Au+Au |l| and S+S at RHIC and LHC. 



Au + Au 



RHIC 



LHC 



e % (GeV/fm 3 ) 
T (GeV) 
to (fm) 



61.40 
0.668 
0.25 
0.34 
0.068 



425 
1.02 
0.25 
0.43 
0.086 



S + S 



RHIC 



LHC 



u (GeV/fm 3 ) 
T (GeV) 
to (fm) 



24.3 
0.531 
0.25 
0.34 
0.068 



170 
0.811 
0.25 
0.43 
0.086 
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FIG. 1. Survival probability of directly produced J/ip at 
RHIC energies due to screening of colour interaction (solid 
curve) and gluonic dissociation in quark gluon plasma. The 
dashed curve gives the latter with inclusion of formation time 
of the charmonium while the dot-dashed curve gives the same 
with the assumption that a fully formed J/tjj is available at 
t — Ti when the plasma is formed. 
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FIG. 3. Same as Fig. 1 for T at LHC energies. The Debye 
screening is absent for the initial conditions fl2| used here. 
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FIG. 2. Same as Fig. 1 at LHC energies. 



